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ABSTRACT 

We present a semi-analytical approach for modeling circumstellar emission from ro- 
tating hot stars with a strong dipole magnetic field tilted at an arbitrary angle to 
the rotation axis. By assuming the rigid- field limit in which material driven (e.g., in 
a wind outflow) from the star is forced to remain in strict rigid-body co-rotation, we 
are able to solve for the effective centrifugal-plus-gravitational potential along each 
field line, and thereby identify the location of potential minima where material is 
prone to accumulate. Applying basic scalings for the surface mass flux of a radiatively 
driven stellar wind, we calculate the circumstellar density distribution that obtains 
once ejected plasma settles into hydrostatic stratification along field lines. The result- 
ing accumulation surface resembles a rigidly rotating, warped disk, tilted such that its 
average surface normal lies between the rotation and magnetic axes. Using a simple 
model of the plasma emissivity, we calculate time-resolved synthetic line spectra for 
the disk. Initial comparisons show an encouraging level of correspondence with the 
observed rotational phase variations of Balmer-line emission profiles from magnetic 
Bp stars like sigma Ori E. 

Key words: stars: magnetic fields - stars: rotation - stars: mass-loss - stars: emission- 
line - stars: chemically peculiar - stars: early-type 
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1 INTRODUCTION 

High resolution images of the solar corona provide vivid evi- 
dence of how the complex solar magnetic field can structur e 
and confine coronal plasma (e.g. jDel Zanna fc Masonl200ffi . 
In other cool, solar-type stars similar complex magnetic 
structuring of a hot corona is inferred indirectly through ro- 
tational modulation of the underlying chromospheric emis- 
sion network, and by year-to-decade timescale modulations 
thoug ht to be analogues of the solar mag netic activity cycle 
(e.g.. IWilsonlll9 78: Baliunas ct al. 1995). By contrast, in a 
subset of hotter, early-type (O, B, and A) stars, spectropo- 
larimetric measurements provide quite direct evidence for 
relatively strong, stable, large scale magnetic fields, of order 
1-10 kG, and generally characterized as a dipole with some 
arbitrary tilt relative to the rotation axis. Instead of a hot 
corona, such stars often exhibit hydrogen Balmer emission 
associated with relatively cool material at temperatures of 
ca. 20, 000 K, comparable to the stellar effective tempera- 
ture. The present paper develops a Rigidly Rotating Mag- 
netosphere (RRM) model for this emission, based on the 
notion that material in the star's radiatively driven stellar 



wind is channeled and confined into co-rotating, circumstel- 
lar clouds by a strong, rigidly rotating dipole field. 

Following the pioneering det ection of stron g fields in 
the chemically peculiar Ap stars (Babcock 1958)), observa- 
tions in the mid- and late-1970s revealed similar magnetic 
fields in both the late B-typ e helium- weak stars (e.g., a 
Cen - IWolff fc Morrison! Il97i) . and th e earlier (types B0- 
B2) h elium-strong stars (e.g., a Ori E - lLandstreet fc Borral 
Il97ct) . More recently, more moderate magnetic fields have 
uission-line stars (e.g., /3 Cep - 
), slowly-pulsating B stars (e.g., £ Cas 
and O-type stars (e.g., 1 Ori C - 



iHenrichs et al 


2000) 


- iNeiner et al 


I2003T 



Many of the Bp stars 1 exhibi t both spectroscopic 
and photometric varia bility (see, e.g.. |Pedersen fc Thomsenl 
Il977t iPedersenl Il979fl . strongly correlated with changes 
in circular polarization arising from their magnetic fields. 
These variations have been interpreted in terms of the 
same 'oblique rotator' conceptual framework that is applied 
to the Ap stars: atmospheric stabilization via a magnetic 
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1 By which we refer to both the helium-weak and the helium- 
strong, chemically-peculiar B-type stars. 
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field allows elemental diffusion to generate surface abun- 
dance anomalies, whose axis of symmetry is parallel to the 
magnetic axis, and t herefore inclined to the rotation axis 
jMichaud et alJll98ll) . 

Often, the variability seen in Bp stars manifests itself 
in circumstellar as well as photospheric diagnostics. Per- 
haps the best-studied example is the B2p helium-strong star 
cr Ori E, which shows Ha shell-like emission varying on the 
same 1.19 dy rotation period as pho tospheric absorption pro- 
files and photometric indices (see iGroote fc Hungerl Il982t 
and references therein) . From studies of this star, and from 
inves tigations of other Bp stars that exhibit sim ilar emission 
(e.g., IShore fc Brownl Il99(j IShore et alJll99Ch . a common 
observational picture has emerged of circumstellar plasma 
confined into tori or clo uds, and for ced into co-rotation, by 
a strong magnetic field JShorelll993l) . 

In the case of a Ori E, the material responsible for both 
the variable Balmer emission and the ec lipse-like behaviour 
seen in photometric light curves fe.g.. iHesser et allll97fif l 
appears to be concentrated at the intersection between the 
rotational and magnetic equators iGroote fc Hungerlll982l: 
iBolton et, aDll987l : Tshort fc Boltonlll994h . An obvious can- 
didate for imposing such structure is the centrifugal accel- 
eration arising from magnetically enforced co-rotation; not 
only can this force lead to the required breaking of symmetry 
about the magnetic axis, it can also furnish the outward lift 
necessary for confining plasma toward the tops of magnetic 
loops. 

This overall scenario is somewhat related to the Mag- 
netically Confined Wind Sh ock (MCWS) model proposed by 
iBabel fc Montmerld dl997allbh to explain the X-ray emission 
from the AOp star IQ Aur and the 07pe star 9 1 Ori C. How- 
ever, their model focuses on the wind collision shocks that 
can produce hot, X-ray emitting gas at the top of closed 
loops, and does not follow the fate of the radiatively cooled, 
post-shoc k material. Magnetohydrody nami cal (MHD ) simu- 
lations bv lud-Doula fc OwockifcOOa) and lud-Douial (120031) 
indicate that, without any rotational support, this material 
simply falls back along the field line to the loop footpoint. 
Nonetheless, recent MHD simulations of the MCWS scenario 
applied to 8 1 Ori C have been quite successful in r eproducing 
its observed X-ray properties IjOagne et alJ 2004). 

Unfortunately, a similar MHD simulation test of the 
magnetocentrifugal confinement scenario is much more dif- 
ficult to carry out. The strong magnetic fields characteristic 
of Bp stars, coupled with the relatively low densities asso- 
ciated with their lower mass loss rates, imply a very high 
Alfven speed. As a result, the time step required to ensure 
numerical stability, via the Courant-Friedrichs-Lewy crite- 
rion, becomes quite short. This makes it very expensive to 
calculate an MHD model spanning the timescales (~ days) 
of interest, even for the relatively simple, two-dimensional 
axisymmetric case of a dipole aligned with the rotation axis. 
For the more general, tilted-dipole case that would apply to 
a Ori E and other magnetic hot stars, the three-dimensional 
nature of the system makes full MHD simulation impracti- 
cal. 

However, in the strong-field limit, an alternative ap- 
proach becomes viable. Under the assumption that the field 
is sufficiently strong so as to remain completely rigid, the 
plasma moves along trajectories that are prescribed a priori 
by field lines that co-rotate with the star. This reduces the 



overall three-dimensional modeling of circumstellar material 
into a series of one-dimensional problems for flow evolving 
under the influence of an effective gravito-centrifugal poten- 
tial. 

iMichel fc Sturrockl lll974ll used such an approach to 
model the magnetosphere of Jupiter, arguing that exo- 
spheric material tends to accumulate in minima of the ef- 
fective potential, occurring along field lin es that pass near 
and above the geostationary orbital radius. Nakaiima 
demonstrated how the same approach can be applied to 
the circumstellar mater ial of oblique rotat or stars such as 
a Ori E. More recently. IPreuss et alJ l|2004h have presented 
an alternative formulation of the strong-field limit, using the 
condition of force balance tangential to field lines to map 
out the complex surfaces on which circumstellar material 
can accumulate. Due to the interrelation between force and 
potential, the latter treatment is entirely equivalent to the 
prior studies based on effective potential minimization. 

In the present paper, we use these studies as the foun- 
dation on which we build the RRM model. Sec. [5] conducts 
a detailed review of the effective-potential formulation for 
the strong-field limit; this review serves both to establish a 
more- rigor ous footing for th e analyses bv lMichel fc Sturrockl 
l|l974h and lNakai imal 11285), and as a basis for the develop- 
ments presented in subsequent sections. Using this formula- 
tion, we examine how the loci of effective-potential minima 
define a likely accumulation surface for circumstellar mate- 
rial (Sec. [jJJ. We then extend our analysis to a full RRM 
model for the circumstellar material, including its hydro- 
static stratification around the potential minima (Sec. |1J, 
its build-up by feeding from the star's wind outflow (Sec. 0, 
and its associated circumstellar line emission (Sec.|SJ. The 
main text concludes (Sec.0 with a comparison of our anal- 
yses with those from previous studies, and with a brief sum- 
mary of results (Sec.|HJ. Finally, the Appendix provides sup- 
porting analyses of the ultimate centrifugal breakout of ac- 
cumulated material against the limited confining effect of a 
finite-strength magnetic field. 



2 THE STRONG-FIELD LIMIT 
2.1 Basic Principles 

In developing a model for the strong-field limit, we adopt two 
basic assumptions. The first is the 'frozen flux' condition of 
ideal MHD, in which plasma is constrained to move along 
magnetic field lines. The second is that these field lines are 
both rigid and time-invariant in the frame of reference that 
rotates at the same angular velocity £1 as the star. Together, 
these assumptions lead to a picture of plasma moving along 
trajectories that are fixed in the co-rotating frame, these 
trajectories being none other than the guiding magnetic field 
lines. 

To develop an understanding of how plasma is chan- 
neled along the rigid field lines, let us first consider the case 
of a solitary parcel launched ballistically along one such line. 
The total instantaneous vector acceleration g tot experienced 
by this parcel, in the co-rotating frame, may be broken down 
into separate components, viz. 

Stot =9 + 9ma g - 2ft x r - ft x (O x r); (1) 
here, g is the acceleration due to gravity, and g mag is that 
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due to the magnetic Lorentz force. The last two terms in 
this expression arise from the inertial Coriolis and centrifu- 
gal forces due to the rotation of the reference frame; fl is the 
vector angular velocity describing this rotation, with mag- 
nitude \tt\ — Q, while r is the position vector of the plasma 
parcel, with its time derivative r giving the corresponding 
velocity vector. 

At any time, the location of the parcel on its respective 
field line may be specified by the arc-length distance s from 
some arbitrary fiducial point. The temporal evolution of this 
field line coordinate is governed by the equation 

d 2 s _ .. „ 

dp = » = fftat ' et. (2) 
where 



(3) 



is the unit vector tangent to the parcel's trajectory. By our 
assumptions above, this trajectory is always directed along 
the field line; therefore, this vector is given by 



et 



B 

~B' 



(4) 



where B is the local magnetic field vector, of magnitude 
B = \B\. 

The equation of motion indicates that the dynam- 
ics of the parcel are governed solely by the component of 
g tot directed along its trajectory. The components of g tot 
perpendicular to e t have no effect on these dynamics: while 
they supply the centripetal acceleration necessary to change 
the direction of the parcel's space velocity r, they leave its 
speed s unchanged. Since both the Coriolis acceleration and 
the Lorentz acceleration 

<? mag = ^(V x B) x B (5) 

are perpendicular to e t (see eqns. [3]and2J, it follows that 
the equation of motion does not depend on the appearance 
of these terms in the expression Q for g tot ; accordingly, we 
find that 



s = [g — fl X (f2 X r)] ■ et. 



(6) 



The gravitational and centrifugal terms in the brackets arise 
from conservative forces; therefore, they may be expressed 
in terms of the gradient of an effective potential "!>, such that 
the equation of motion becomes 



s = ~(V$) ■ e t 



(7) 



(8) 



Recognizing the right-hand-side as the directional derivative 
of $ along the field line, we have 

s = -dT 

This result is very instructive: it tells us that although the 
plasma parcel follows a three-dimensional curve r = r(s) 
through space, its motion is governed by a potential func- 
tion 3>(s) that arises from sampling $(r) along this curve. 
Throughout, we term this single-variable function the field 
line potential. 

If the field line potential $(s) exhibits an extremum, so 

that 



as 



at some point, then the plasma parcel can remain at rest, 
with no net forces acting upon it; the components of the 
gravitational and centrifugal forces tangential to the field 
line are equal and opposite, while those perpendicular to the 
field line are balanced by the magnetic tension. Whether the 
parcel can remain at such an equilibrium point over signifi- 
cant timescales (viz., multiple rotation periods) depends on 
the nature of the extremum. In the case of a local maximum, 
where d 2 $/ds 2 = $" < 0, the equilibrium is unstable, and 
small displacements away from the extremal point grow in 
a secular manner. This is what lud-Doula fc OwockTll2002h 
found in their MHD simulations of wind outflow from a non- 
rotating star, in the case of a moderately strong magnetic 
field (their n, = 10). Since the effective potential in the ab- 
sence of rotation is just the gravitational potential, the tops 
of magnetic loops are local maxima of the field line poten- 
tial. Therefore, although plasma at loop tops is supported 
against gravity by magnetic tension, it is unstable against 
small perturbations, and - as the MHD simulations show - 
it eventually slides down one or the other side of the loop 
toward the stellar surface. 

In the converse situation, where the extremum in the 
field line potential is a minimum, with $" > 0, the equi- 
librium is stable: any small displacement along the local 
magnetic field line produces a restoring force directed back 
toward the equilibrium point. Such minima represent ideal 
locations for circumstellar plasma to accumulate. Because 
these potential minima can occur on more than a single 
field line, the accumulation is not at an isolated point in 
space, but rather is spread across one or more surfaces de- 
fined by the loci at which both $' = and 4>" > 0. Material 
that collects on these accumulation surfaces forms a mag- 
netosphere that is at rest in the co-rotating reference frame; 
when viewed from an inertial frame of reference, this mag- 
netosphere appears to rotate rigidly with the star, hence the 
name chosen for the RRM model. 

As a brief aside, it is readily demonstrated that the fore- 
going potential-based analysis is entir ely equivalent to th e 
force-based formulation presented by IPreuss et alJ i2004l) . 
For instance, the condition $' = for an equilibrium point 
(stable or unstable) may be expressed in the form 



[g - n x (O x r)] • B = 0, 



(10) 



which comes from combining eqns. (HIHIIl with eqn. <QJ. This 
expression can be recognized as the exact same condition 
of force equi l ibrium tangential to the local field line that 
IPreuss et alJ (|2004T l impose in their eqn. (2). 

We turn now to examining the effective potential 4>(r), 
which determines the potential 5?(s) along each field line. 
Within the Roche limit, where the star is assumed to be so 
centrally condensed that it may be treated as a point mass, 
this effective potential is given by 



$( r ) = 



GM» 



l n 2 2 • 2 

Si r sin i 

r 2 



(11) 



(9) 



here, G is the gravitational constant, M* the stellar mass, 
and r and 6 are the radial and colatitude coordinates cor- 
responding to the position vector r, in the spherical-polar 
system aligned with the rotation axis. Let us introduce the 
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Kepler co-rotation radius 2 



*(x,z) [j3 = 0°] 



V fi 2 



(12) 



at which the gravitational and centrifugal forces balance in 
the equatorial plane; then $ may be written as 



<f>(r) = 



GM, 



1 lj . . 2 



(13) 



where £ = r /rjc is the radial coordinate in units of the Kepler 
radius. This latter form is convenient, because the minima 
of $(s), along each field line, occur in the same location 
as the corresponding minima defined by the dimensionless 
potential 



*(r) 



GM, 



Mr) 



1*2 ■ 2 a 

— t sm 0. 

2 s 



(14) 



The advantage of working with this dimensionless poten- 
tial is that i t is i ndependent of the rotation rate Q; as 
iPreuss et alJ <l2004t) demonstrate, a similar conclusion can 
be reached in the force-based formulation of the problem. 
Accordingly, for each magnetic field configuration, we only 
need solve once for the accumulation surfaces where plasma 
can remain at rest in stable equilibrium; this solution can 
then be mapped onto a specific rotation rate by transform- 
ing the radial coordinate from £ back to r. 

Looking at the form of the dimensionless potential ^(r) 
introduced above, we can identify two regimes. When r is 
much smaller than the Kepler radius rx, such that ^ < 1, 
this potential is spherically symmetric about the origin, and 
increases outwards. Conversely, when the distance from the 
rotation axis greatly exceeds the Kepler radius, such that 
£sin# ^> 1, the potential is cylindrically symmetric about 
the same axis, and decreases outwards. As we demonstrate 
in the following sections, it is in this second regime that the 
field line potential <3?(s), and its dimensionless equivalent 
^(s), exhibit the minima near which circumstellar plasma 
accumulates. 



2.2 The aligned dipole configuration 

In the foregoing discussion, we argue that circumstellar 
plasma accumulates on the surfaces defined by minima of 
the field line potential. We now consider the geometry of 
one such accumulation surface in the simplest of all con- 
figurations, that of a centred dipole magnetic field aligned 
with the star's rotation axis. Defining a Cartesian coordinate 
system (x, y, z) with origin at the star's centre and rotation 
taken along the z-axis, Fig.^shows an x vs. z contour map 
of the dimensionless effective potential $(r) (see eqn. I IH 
in the y = plane. Superimposed over the map are four 
curves, each following the parametric equation for a dipole 
field line, 



r ■ 2 

£ = 7 sin 



(15) 



Here, the parameter 7 specifies the summit radius (in units 
of nc) of the field line; in the present case of Fig. the 
7 = y/2/3, y/2, 2 and 4 lines are plotted. The significance of 

2 Applied to Earth, the Kepler radius is equivalent to the orbital 
radius of a geostationary satellite. 
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Figure 1. A contour map of the dimensionless effective potential 
over the y = plane; white regions corresponds to > —0.6, 
black regions to < —3, and the intermediate gray levels are 
spaced in increments = 0.15. Drawn over the map are four 
selected field lines for a magnetic dipole aligned with the rota- 
tion axis; these lines have summit radii 7 = y2/3 (solid), \/2 
(dashed), 2 (dotted) and 4 (dot-dashed), and are oriented at az- 
imuthal angles <f> = 0° (right) and <j> = 180° (left). Shown beneath 
is the dimensionless potential along each of the = 0° lines, 
plotted as a function of 9; note that the scale of the ordinate is 
logarithmic. 



the first value is discussed below. Beneath the contour map, 
we plot the dimensionless potential 



no) 



1 2 • 6 , 

-7 sm ( 



(16) 



7sin 2 2 

this being $(r) sampled along the dipole trajectory 1151 . 
for each of the four field lines. For simplicity, we chose the 
colatitude 8 as the independent variable in the above expres- 
sion, rather than the usual field line coordinate s. However, 
noting that the two are related via the differential equation 



ds 
AO 



■■ r K 7 



sin 6y/l + 3 cos 2 9, 



(17) 



it is clear that, everywhere away from the poles, s varies 
monotonically with 9. 

Inspecting the *&(#) data for the three outer field lines 
(7 = \/2, 2, 4) shown in Fig. □ minima can be seen at the 
stellar equator (9 = 90°). Since the aligned dipole configura- 
tion is symmetric about the z axis, we can conclude that the 
accumulation surface takes the form of an equatorial disk, 
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with its normal pointing along the rotation axis. However, 
because a potential minimum does not occur along the in- 
nermost field line, it is evident that this disk does not extend 
to the origin, but instead must terminate at some finite ra- 
dius. To determine this inner truncation radius, we observe 
that 



d 2 # d9 d 2 6> d\g 
1(P + ds ds2 d# 



(18) 



for 9 — 90°, where in the second line we make use of the 
identity £ = 7 within the equatorial plane (cf. ean. 1151 . We 
recall that 9" must be positive in order for an extremum 
(ty 1 = 0) to constitute part of an accumulation surface; ac- 
cordingly, the inner truncation radius is given by 



6 



0.87 



(19) 



at which changes from being positive (£ > £i) to negative 

« < 6). 

The reason for our choice of 7 = \/2/3 = £s for the in- 
nermost field line in Fig.^should now become apparent: it 
ensures that this particular line exactly intersects the inner 
edge of the accumulation disk at £ = £i. Plasma at the sum- 
mit of this line is therefore in neutral equilibrium, whereby 
small displacements away from the equator produce no net 
force (to first order in the displacement) along the field line, 
either away from the equilibrium point or toward it. This is 
evident in the lower panel of Fig. from the flatness of this 
field line's dimensionless potential at 9 — 90°. 

Throughout the region in the equatorial plane between 
the truncation radius £ = £j and the Kepler co-rotation ra- 
dius 3 £ = 1, magnetic tension supports accumulated mate- 
rial against the net inward pull caused by gravity exceeding 
the centrifugal force. Beyond this region, when £ > 1, the 
centrifugal force surpasses gravity, and the effect of magnetic 
tension then becomes to hold material down against the net 
outward pull. 

Clearly, the interplay between gravitational and cen- 
trifugal forces has a different significance in the RRM case 
than it does for a Keplerian disk; in the latter, material 
at each radius orbits the star at a velocity whereby both 
forces are in exact balance. Such complete force balance is 
not required in an RRM, inasmuch that magnetic tension 
can a bsorb any net resu ltant force perpendicular to field 
lines dPreuss et al.l l2004'l . Only the tangential components 
of the forces are required to be in balance, so as to produce 
an equilibrium that is stable against small displacements 
along fie ld lines - a point r e cogni zed in the original treat- 
ment bv lMichel fc Sturrockl (Il974f) . although these authors 
employed a more-geometrical approach to arrive at the same 
conclusion. 



2.3 The tilted dipole configuration 

Up until now, we have dealt with the trivial case of a dipole 
field aligned with the rotation axis. However, the re is noth- 
ing that restricts us to such simple systems; as iNakaiimal 



(1985) first demonstrated, the rigid-field approach we have 
presented can be applied to arbitrary magnetic configura- 
tions, so long as the effective potential along each field line 
can be computed and minimized. 

In the present section we now consider the oblique ro- 
tator configuration, where a dipole is inclined at an angle j3 
to the rotation axis. For such a geometry, eqn. 111! still de- 
scribes the effective potential in the co-rotating frame, but 
the field lines of the tilted dipole now follow the parametric 
equation 



£ = 7 sin 6, 



(20) 



where 9 is the colatitude coordinate in the frame of refer- 
ence aligned with the magnetic axis. To relate this magnetic 
reference frame back to the rotational one, we adopt the 
convention that the former is obtained from the latter by 
rotating by an angle (3 about the Cartesian y axis 4 . With 
this convention, colatitudes in the two reference frames are 
related to one another via 



sin 2 9 = sin 2 9 sin 2 6 + 



sin (3 cos 9 + cos (3 sin 9 cos (f> ) , (21) 



were <j> denotes the azimuthal coordinate in the magnetic 
frame. 

The latter expression may be used to eliminate the 
sin 2 9 term from eqn. (1141 . allowing us to express the di- 
mensionless effective potential along each field line as 



*(*) = 



7sm 



1 2 ■ 4 2 \ ■ 2 o -27, 

-7 sm a sin 9 sin <p + 



sin f3 cos 9 + cos /3 sin 9 cos < 



(22) 



here, we have also used eqn. (1201 to eliminate £. It is straight- 
forward to derive an expression for the derivative VP' of this 
field line potential; however, in contrast to the aligned dipole 
configuration, the equation ^' — for the extrema of this 
potential no longer admits algebraic solutions. But in the x- 
z plane that contains both the magnetic and rotation axes, 
we can still illustrate these minima using the same graphical 
approach as in Sec. 12.21 

Figure|5|shows similar effective potential plots to Fig.0 
but for dipole field configurations tilted at angles (3 = 30°, 
60° and 90° to the rotation axis. Focusing initially on the 
first two cases, we note that a single potential minimum 
occurs along the outer three field lines when f3 — 30°, and 
along the outer two for (3 = 60° . In each case, the minima are 
situated at approximately the same colatitude, which falls 
somewhere between the magnetic and rotational equators: 
9 » 80° for 13 = 30°, and 9 » 70° for (3 = 60°. This bisection 
of the equators arises because of competition between the 
two misaligned symmetry axes, magnetic vs. rotational. 

Looking now at the (3 = 90° case, we can see that two 
minima - albeit shallow ones - occur in the potential along 
the outermost field line, at equal distances above (9 « 55°) 
and below (9 ~ 125°) the magnetic equator. This shows 
that the accumulation surfaces of tilted configurations can 



3 Apparent in the figure as the twin saddle points at (x, z) = 
(±1,0) r K . 



4 We assume that /3 > corresponds to a clockwise rotation, 
when looking out from the origin along the positive y axis. 
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Figure 2. As with Fig.Q except that a tilted dipole field is assumed, at angles /3 = 30°, 60° and 90° to the rotation axis. 



be significantly more complex than the simple disk found 
for the aligned dipole; indeed, as we demonstrate in the fol- 
lowing section, there is no guarantee even that the surfaces 
are made from a single contiguous sheet. iNakaiimal (I1985T) 
ov erlooked such situat ions, and it was not until the study 
bv lPreuss et alJ l)2004) that this possibility became known. 

To conclude this section, we draw attention to a limi- 
tation of the graphical approach we have used to illustrate 
the accumulation surfaces: although these surfaces are in- 
herently three dimensional, our approach is restricted to 
plotting the effective potential and magnetic field lines over 
a two dimensional slice through the system. This is not a 
problem for the aligned dipole shown in Fig. Q since rota- 
tional symmetry ensures that all slices containing the polar 
axis are identical. However, this symmetry is absent from the 
tilted field configurations plotted in Fig. El meaning that the 
figure cannot indicate the nature of the accumulation sur- 
faces outside of the x—z plane that contains the misaligned 
rotation and magnetic axes. 



3 ACCUMULATION SURFACES 

To overcome these limitations of two-dimensional plots, 
let us now use perspective images to show the full three- 
dimensional form of the accumulation surfaces. For the same 
aligned and tilted dipole magnetic field configurations intro- 
duced above, Fig. ^illustrates these as surfaces illuminated 
by an artificial parallel light source from the observation 



point. The vertical, single-headed arrows denote the rotation 
axis, while the double-headed arrows shown in projection at 
various orientations represent the magnetic axis. For each of 
the four values of the tilt angle (3 = {0°, 30° , 60°, 90°}, the 
accumulation surface is shown from three different observa- 
tion points, situated at the same inclination i = 60° to the 
rotation axis, but having differing azimuthal angles (f> — 0° , 
45° and 90°. Although the full surfaces formally extend to 
arbitrarily large radii, for the illustration the outer edge is 
truncated by omitting regions threaded by field lines with 
summit radii greater than 7 = 12. 

For the aligned field (/3 = 0°) case, the accumulation 
surface is a simple disk lying in the plane of the rotational 
and magnetic equators, and so appears identical from all az- 
imuths. The hole at the centre reflects the lack of potential 
minima inside the inner truncation radius (Sec. 12.21 . For 
the /3 = 30° case the surface is titled, with a mean normal 
vector between the two symmetry axes, in a direction con- 
sistent with the 6 w 80° angle found in Sec. 12.31 Although 
not obvious in Fig. |3 the surface is not strictly planar, but 
has a slight warp. 

For the greater tilt of the j3 — 60° configuration, this 
warping becomes more apparent. While still shaped approx- 
imately like a disk, the regions nearest the intersection with 
the plane formed by the two axes are warped away from the 
rotation axis. Physically, this arises because the centrifugal 
force vanishes along the rotation axis; this force being crucial 
to forming the potential minima that make up the accumu- 
lation surface, it follows that there is a 'zone of avoidance' 



,3 = 0°, O = O° 
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,8 = 30°, O =O° 




Figure 3. Accumulation surfaces for aligned (/3 = 0°; left-hand column) and tilted ({3 = 30°; right-hand column) dipole magnetic fields; 
the surfaces are viewed at an inclination i = 60° to the rotation axis, and from azimuths (top-to-bottom) <j> = 0°, 45° and 90°. In each 
plot, the rotation and magnetic axes are shown as single- and double-headed arrows, respectively; to indicate the scale of the plots, a 
bar with a length of one Kepler radius (rx) is shown in the bottom right-hand corner of each. 
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Figure 3 - continued As before, except that the tilted configurations for (i 
are shown. 



= 60° (left-hand column) and /3 = 90° (right-hand column) 
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around the rotation axis, inside of which the surface can- 
not exist. An additional, remarkable feature of the (3 = 60° 
case is the appearance of a pair of secondary accumulation 
surfaces, situated in each hemisphere between the magnetic 
and rotation axes. These secondary surfaces, which we term 
leaves, are a consequence of the appearance of an additional 
minimum in the effective potential along a particular bundle 
of magnetic field lines. 

From the analysis in Sec. 12.31 it might appear that such 
two-minima scenarios are restricted to the perpendicular 
(/3 = 90°) configuration; however, the appearance of the 
le aves in this intermed iate-tilt case proves otherwise. In fact, 
as lPreuss et alJ J2004T) have demonstrated, leaves 5 occur in 
all configurations other than the aligned field one, but are 
situated at larger and larger radii as j3 decreases toward zero. 
In the present case with j3 = 60°, the leaves are threaded by 
magnetic field lines for which 7 > 10; this explains why the 
middle panel of Fig. [5] which plots field lines up to 7 = 4, 
does not exhibit the second potential minimum associated 
with the appearance of a leaf. 

Turning finally to the perpendicular field (/3 = 90°) con- 
figuration, we see the ultimate product of the disk warping 
and leaf formation. The accumulation surface now takes the 
form of a partial disk lying in the magnetic equator, inter- 
sected by an opposing pair of truncated cones aligned with 
the rotation axis. These cones have an opening half-angle of 
tan" 1 2~ 1,/2 « 35.3° at asymptotically-large radii, a value 
that can be derived by setting the first derivative of 9(9) 
(cf. cqn. 1221 to z ero, a nd then solving for as 7 — > 00 (see 
also iPreuss et aJll2004 . To understand the unusual geome- 
try of the perpendicular configuration, note that if j3 were 
to depart slightly from 90°, then the half of each cone that 
lies between the two axes would split off from the main ac- 
cumulation surface, and take the form of a leaf. Therefore, 
we can recognize the cones as being formed from a merger 
between the warped disk and the leaves. 

The geometrical complexity of accumulation surfaces, 
even in t he relatively-simp le case of a tilted dipole, was un- 
known to lNakaiimal jl985T) : he had to rely on simple particle- 
based maps for visualization (see his figs. 2 & 3), and was 
unaware of the possibility of leaves or of the truncated-cone 
configuration occurring at /3 = 90°. Only in recent years 
have computers become sufficiently powerful that the visu- 
alization of the surfaces is a relatively straightforward proce- 
dure. However, there still remains the question of how closely 
a physical system would resemble an accumulation surface. 
The answer depends on the nature and distribution of the 
matter that populates the surfaces and renders them visible 
or detectable. While the accumulation s urfaces presented in 
this section, and bv lPreuss et alJ 1120041) . furnish a geomet- 
rical picture of where circumstellar material can accumu- 
late, they provide no indication of how much material does 
accumulate, nor of its physical conditions - density, tem- 
perature, opacity, emissivity, etc.. We address these issues 
in the following sections, where we derive the distribution 



5 Or 'stable chimney regions' in the parlance of IPreuss et al.| 
J2004t) . the chimney being the rotation-axis aligned surface de- 
fined by $' = 0, that is composed of both stable and unstable 
equilibrium loci. 



of circumstellar gas, and then use this to calculate the line 
emission from an RRM model. 



4 HYDROSTATIC STRATIFICATION 

The various processes that could fill effective potential wells 
with material may generally be quite dynamic and variable; 
but over time, it seems likely that most such material should 
eventually settle into a nearly steady, static state. Along any 
given field line that intersects with one or more of the poten- 
tial minima that define accumulation surfaces, the relative 
distribution of material density p is then set by the require- 
ment of hydrostatic stratification within the field line poten- 
tial, 



dp 



d$ 



ds ^ ds 

where the gas pressure is given by the ideal gas law, 
pkT 
I* 



V 



(23) 



(24) 



Taking for simplicity the temperature T and mean molec- 
ular weight p, to be constant, we can solve for the density 
distribution along the field line as 



p m e 



-/i[*(s)-* m ]/kT 



(25) 



where p m = p(s m ) and <& m = "I?(s m ), with s m the field line 
coordinate at the potential minimum 6 . In our application to 
magnetic hot stars, we assume radiative equilibrium should 
set the circumstellar temperature to be near the stellar ef- 
fective temperature. This implies a thermal energy per unit 
mass kT/p that is much smaller than the variation in the 
corresponding potential energy $, which by eqn. 11131 is typ- 
ically comparable or greater than the gravitational binding 
energy at the Kepler radius. By eqn. I25L we thus expect 
that most of the material should be confined to a relative 
narrow layer near the potential minimum at s = s m . 

A Taylor series expansion of the potential about this 
minimum gives 



*(s) = $m + 



1 - 



(26) 



where we have used the fact that, by definition, "1?^ = 0. Ac- 
cordingly, in the neighborhood of the minimum, the density 
distribution 125H may be well approximated by 



PmS 



where the RRM scale height is 



2kT/p 



L in 



2kT/p 
GM,/r K 



1 



(27) 



(28) 



In the latter equality, the first square root is of the ratio 
between the thermal energy and the gravitational binding 
energy at the Kepler radius, while the second root gives 
the curvature length of the effective potential. Overall, we 
thus see that the effect of a finite gas pressure is to support 
material in a nearly Gaussian stratification on either side of 
an accumulation surface. 



6 Here and throughout we use the subscript "m" to denote quan- 
tities evaluated on the accumulation surface s = s m . 
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A similar Gaussian stratification is found in models for 
Keplerian disks; however, the scale height in that case grows 
as th e three-halves power of the d istance from disk centre 
(e.g., iHummel fc Hanuschiklll997fl . leading to a disk that 
flares outward with increasing radius. In contrast, the RRM 
scale height approaches a constant value far from the ori- 
gin; for example, in the case of the aligned dipole configura- 
tion, eqn. 118H gives = 3/rx when £ 3 3> 1, yielding the 
asymptotic scale height 



I'K 



2fcT//j 
3GM,/r K 



(r > r K ) 



(29) 



(compare with eqn. 5 of lNakaiimalll985h . Since the ratio of 
thermal energy to gravitational binding energy is typically 
very small, we thus see that the disk is indeed geometrically 
thin, with a scale height fe m that is much smaller than the 
characteristic Kepler radius tk- 

A convenient way to characterize the properties of such 
a thin accumulation layer is in terms of its surface density. 
For field line with a projection cosine fim to the surface nor- 
mal, the associated local surface density <j m can be obtained 
by integration of eqn. I|27[l over the Gaussian hydrostatic 
stratification, 



V^PmUmhn 



(30) 



For the application here to magnetic hot stars, we next de- 
velop a specific model for the global distribution of this sur- 
face density as being proportional to the rate of material 
build-up by the loading from the star's radiatively driven 
wind. 



5 MASS LOADING OF ACCUMULATION 
SURFACES 

5.1 The Accumulation Rate for Surface Density 

The high luminosity of hot stars is understood to give rise 
to a radiatively driven stellar wind; in a magnetic hot star, 
this wind provides a key mechanism to load mass into the 
effective gravito-centrifugal potential wells around the accu- 
mulation surfaces. 

For a radiatively driven wind in t he presence of a mag- 
netic field. lOwocki fc ud-Doulal (I2004T) derive an expression 
for the mass flux density at the stellar surface, in terms of 
the spherical mass loss rate M predic ted by the standard 
CAK wind model (ICastor et al.lll975l) . For a dipole flux- 
tube bundle intersecting the stellar surface r = _R* with a 
projection cosine /t», and having a cross-sectional area dA t , 
the rate of mass increase is 



2/uM 



dA* 



(31) 



where the factor 2 accounts for the mass injection from two 
distinct footpoints. 

In a highly supersonic stellar wind, the collision of mate- 
rial from opposite footpoints leads to strong shocks that heat 
the plasma initially to te mperatures of millions of deg rees; 
this dogma, advanced bv lBabel fc Montmerlel l|l997al lb]) in 
their MCWS paradigm , has been amply confirme d both by 
MHD modeling (e.g.. lud-Doula fc Owocki 120021) . and by 



analysis of the observed X-ra y emission from t he super- 
heated post-shock plasma (e.g.. lGaene et al Eventu- 
ally, the plasma cools radiatively back to temperatures near 
the stellar effective temperature of a few times 10 4 K. As 
discussed already in Sec. [I] at such temperatures material 
trapped within the effective potential well will quickly settle 
into a relatively narrow hydrostatic stratification centred on 
the potential minima. 

For the moment, let us consider the relatively simple, 
common case that there is a single minimum at field line co- 
ordinate s m , at which point the flux tube has area dA m and 
intersection cosine /t m with the accumulation surface nor- 
mal. Then the associated rate of increase in surface density 
<r m can be written 



&m = rh/i m /dA m = n„ 



M dA* 



(32) 



2ttRI dA m ' 

From the conservation of magnetic flux, V • B = 0, we have 



dA,B, = dA m B 

m j 

which when applied to eqn. (1321 gives 
M B m 



a m = /i» /in 



2-kRI B, 



(33) 



(34) 



For a dipole field, this declines with radius r as & m oc B oc 

_-3 



5.2 The Time Evolution of the Volume Density 

The above merely gives the rate at which surface density 
increases in the accumulation surface. In the idealized limit 
that the field is arbitrarily strong, the actual surface den- 
sity could thus increase without bound. In reality, for any 
large but finite field, the finite magnetic tension could only 
contain a finite mass of material. As analyzed in the Ap- 
pendix, above some breakout density the net centrifugal and 
gravitational force should overwhelm the tension, leading to 
centrifugal ejection that effectively empties the accumula- 
tion surface. 

This view suggests a simple model in which the local 
surface density builds linearly with the time t since the last 
evacuation, <r m = a m t. Applying such a model to eqns 1301 
and 13411 . we can eliminate the normalizing term p m from 
the equation 12511 for the volume density, to find the time 
evolution of this density as 



p{s,t) 



Mt fj,* 



m - M [*(s)-* m ]/fcT 



(35) 



2tt 3 / 2 Rlhm Bt, 

For field lines that exhibit two potential minima, we divide 
the latter expression by the factor 



r i _(_ [ha 

1 ~ + h~ at 



(36) 



where fv m , Bj, and are the scale height, field strength and 
potential evaluated at the secondary minimum with field line 
coordinate sj n . This factor accounts for the partitioning of 
plasma between the two minima, assuming a free exchange 
of material leads to a common hydrostatic stratification. 

We emphasize that eqn. I|35|l applies only to field lines 
that intersect one or more accumulation surfaces, with one 
or more potential minima. For all other field lines we set the 
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density to zero, reflecting the notion that there can be no 
stable accumulation of material over time. Moreover, even 
for field lines intersecting a surface, we only apply eqn. I|35[l 
up to the bracketing maxima in the potential. Beyond these 
maxima, the decreasing potential suggests an exponentially 
growing density, contrary to the true physical picture of no 
accumulation. We resolve this difficulty by setting the den- 
sity p to zero in the regions beyond the potential maxima 7 . 
This approach creates discontinuities in the densities at the 
maxima themselves, but in most cases ones so small that 
they are unlikely to be of much significance for the overall 
model. 



6 CIRCUMSTELLAR EMISSION 

Let us now apply the expression (1351 for the density dis- 
tribution toward calculating the circumstellar line emission 
from a tilted-dipole RRM configuration. As a typical exam- 
ple case, we adopt a tilt angle /3 = 60°, and set the stellar 
angular velocity fl to 50 percent of the critical rate, 



8GM* 
27R 3 , 



(37) 



at which the surface centrifugal force at the equator would 
balance gravity 8 ; this rotation rate corresponds to a dimen- 
sionless stellar radius = R*/rj< = 0.42. Our choice of pa- 
rameters is loosely g uided by the rotation rate and magnetic 
tilt inferred from the lGroote fc Hunger! il997f) model for the 
surface of a Ori E; however, we make no attempt at fine tun- 
ing, since we are more concerned here with demonstrating 
the capabilities of the RRM model, than with obtaining an 
accurate picture of this particular helium-strong star. 

To specify the temperature T and mean molecular 
weight /i in the model, we introduce the dimensionless quan- 
tity 



kTR* 
fiGM* 



(38) 



which characterizes the ratio of thermal to gravitational 
binding energy at the stellar surface. In the photospheres 
of early-type stars, this ratio is of the order ~ 10~ 3 . Follow- 
ing again the scenario that the circumstellar environment 
remains at a temperature close to photospheric (cf. SecHJ, 
we therefore set e* = 10 -3 throughout. 

We now make the assumption that the plasma volume 
emissivity j\, at a wavelength A, may be characterized by 
the relation 



Ja = j p 2 5(\- Ao); 



(39) 



here, jo and Ao are constants, and 5(. . .) is the Dirac delta 
function. This expression is intended to mirror the process 



7 Unless, of course, a region beyond a potential maximum belongs 
to a neighbouring secondary minimum at s = - in which case, 
the expression for the density remains valid. 

8 Note that, while this expression for Q c is appropriate to a 
centrifugally-distorted star, we have assumed elsewhere, for sim- 
plicity, that the star remains spherical. For the particular choice 
Q = 0.5 n c , the stellar oblateness remains small, with the equato- 
rial surface radius being barely 4 percent larger than that at the 
poles. 



of monochromatic line emission at a rest-frame wavelength 
Ao, arising from the density-squared radiative recombina- 
tion of ionized hydrogen. Integrating the emissivity along a 
given ray passing through the magnetosphere, the observed 
surface intensity of the emission is given by 



T x = 



f°° 2 

jo / p (zo) S(X - A [l + v p (z )/c\) dz a 
Jo 



(40) 



where z a is the distance along the ray from the observer, 
v p is the projection of the local plasma velocity onto the 
ray, and c is the speed of light. For rays intersecting the 
star, this integral must be truncated at the stellar surface, 
to account for the occultation of radiation incident from the 
star's far side. Note that this simple expression does not in- 
clude the emission of radiation by the star itself; therefore, it 
should be taken to represent the notional circumstellar com- 
ponent of the net radiation from the system, with the un- 
derstanding that the corresponding photospheric component 
has been subtracted away. Of course, such an interpretation 
is in itself an approximation, since our emission model does 
not account for episodes when circumstellar plasma tran- 
sits the disk and absorbs stellar radiation over the range 
— v e smi < v p < w e sin i; however, at the level of the present 
investigation, this approximation is entirely sufficient. 

Because the magnetospheric plasma co-rotates rigidly 
with the star, v p may be expressed as 



Vp(z ) = A n e sin i, 



(41) 



where the impact parameter A is the perpendicular distance, 
in units of R, , between the ray and the rotation axis, and v e 
is the equatorial rotation velocity of the star. These quanti- 
ties, and therefore v p too, are independent of z ; hence, the 
intensity may be written in the form 



: jo 5(X - Ao[l + Aw e sini/c]) / p 2 (z )Az . 

Jo 



(42) 



From this latter expression, it can be seen that all plasma 
having the same A - that is, situated on a plane parallel 
to the rotation axis - will radiate monochromatically at the 
same wavelength. 

Applying eqn. 1421 to the known density distribution of 
the j3 — 60° RRM configuration, Fig. 0] shows maps of the 
wavelength-integrated emission intensity 



I\ dA, 



(43) 



extending out to 5 R* in directions parallel and perpendic- 
ular to the projected rotation axis. The observer is situated 
at the same inclination i = 60° that we adopt in Sec. [31 
and at eight differing values of the azimuth <f) , separated by 
uniform increments of 45°. Beneath each map we show the 
corresponding spectrum, in which the spatially-integrated 
emission is plotted as a function of the projected velocity 
Vp. Since we are interested more in the distribution of emis- 
sion than its absolute value, Fig. 0] adopts an arbitrary (al- 
though consistent) normalization for the intensities in both 
the maps and the spectra. 

In Section|3] we demonstrate that the accumulation sur- 
face for a P = 60° tilted dipole takes the form of a warped 
disk. The emission maps in Fig. 01 reveal that the distribu- 
tion of material across this disk is decidedly non-uniform. 
Specifically, the distribution is dominated by two clouds, 
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= 135 c 




v p (v e sin i) 



Figure 4. Maps of the emission from the j3 = 60° RRM model (see text), for eight different values of the observer azimuth O ; darker 
shading indicates greater intensity. In each map, the rotation and magnetic axes are shown as a single- and double-headed arrows, while 
the outline of the central star (whose contributions toward the emission are neglected) is shown by a circle. Beneath each map, we plot 
the corresponding emission spectrum as a function of the projected velocity v p , the latter being measured in units of the projected 
equatorial rotation velocity v e sini of the star. 



located near the inner edge of the disk at the intersection 
between the rotational and magnetic equators. Seen from an 
inertial frame, these clouds appear to rotate synchronously 
with the star; furthermore, their characteristic twin-peaked 
emission spectrum displays temporal variations in the form 
of a double S-wave. 

These findings exhibit an encouraging degree of agree- 
ment with the inferred behaviour of a Ori E. As discussed 
in the introduction, observations of this star indicate that 
circumstellar plasma is concentrated at t he intersection be- 
tween rotational and magnetic equators (IGroote fc Huneeii 
ll982l:lBolton et alJl987tlShort fc Boltonll994) . Without the 
need for any special tuning, beyond the requirement that the 
tilt angle /3 be moderate, the RRM model naturally repro- 
duces such a distribution. It also accounts for the Balmer 



line emission from the star, which - when the measured lon- 
gitudinal field is strongest, corresponding to <f) — 0° and 
<f>o = 180° in Fig.Q]- is observed to exhibit strong peaks sit- 
uated at v p ~ ±3f e sini (Groote, personal communication). 



7 DISCUSSION 

As we discuss in Sec. there have been a number of pre- 
vious studies that have made use of the rigid-field approach 
(cf. Sec. to determining the regions where circumstellar 
material can accumulate (Sec.|SJl. The present paper builds 
on these studies, by presenting a physically-grounded RRM 
model for the steady accumulation of wind plasma in the 
circumstellar environment (Sees. 3] & that is able to 
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make specific predictions regarding the observables associ- 
ated with this plasma (Sec. |SJ . 

Our treatment of the mass loading of accumulation 
surfaces d iffers markedly from the approach advanced by 
INakaiimal lIlQSfiT ) , who fixed the plasma density at each point 
by requiring equal magnetic and kinetic energy densities (see 
his eqn. 13); such a choice is guided by the notion that when 
the density is high enough for the kinetic energy due to 
rotation to dominate the magnetic energy, the field lines 
break open, and any subsequently-added plasma leaks away 
from the system. By contrast, our approach to deriving the 
plasma distribution focuses on the mass accumulation rather 
than on leakage (cf. Sec. 0. While INakaiimal il985T) treats 
the mass leakage as a gradual, quasi-steady process, we view 
it more as an episodic evacuation caused by magnetic break- 
out (see Appendix), which effectively resets the mass accu- 
mulation. 

In actual systems the mass distribution may reflect a 
combination of both perspectives, and there are even other 
alternative frameworks f or treating the problem (see, e.g., 
iMichel fc Sturrocklll974l : lHavnes fc Goertzlll984D . However, 
one particularly favourable aspect of the present model is 
that it can naturally reproduce the plasma concentrations 
at the intersection of rotational and magnetic equators, as 
is infe rred from observa tions of a Ori E. To obtain a similar 
result. INakaiimal il985T) had to make the ad hoc assumption 
that some process of diffusion, across magnetic field lines, 
leads to the redistribution of plasma into the desired config- 
uration. 

We turn now to a brief discussion of the recent paper 
bv lPreuss et alJ <l2004l) . which was published during the fi- 
nal stages of preparation of the present work. These au- 
thors found the same accumulation surfaces as we present 
in Sec. El but using the alternative formulation based around 
consideration of the loci where all forces tangential to field 
lines vanish. They were the first to discover the possibility 
of leaves and the truncated cones occurring at /3 = 90° , and 
even expl ored the case of tilte d dipoles offset from the origin. 
However, iPreuss et alJ (I2004T) stopped short of considering 
the mass loading of the accumulation surfaces, and instead 
focused on their geometrical form. As such, their study does 
not make the kind of specific predictions for observable emis- 
sion line variations that we provide in Sec. |S| 

Another relevant comparison here is to the Magnetically 
Torqu ed Disk (MTD) model proposed by ICassinelli et alJ 
(2002) to explain the circumstellar emission of Be stars. 
Building upon insights from one-dimensional equato- 
rial plane models of magnetically torqued stellar winds 
JWeber fc Davis! Il967t iBelcher fc MacGreeorl Il97d) . this 
analysis centres on an assumed empirical scaling of the az- 
imuthal velocity, which initially increases as a rigid-body law 
out to some peak, and then declines asymptotically with 
angular momentum conservation. When this peak occurs 
above the Kepler radius, the model envisions that the asso- 
ciated torquing of the wind outflow can lead to formation of 
a 'quasi-Keplerian' disk. B ut recent dynamical simulations 
JOwocki fc ud-Doulalliool) indicate that fields marginally 
strong enough to spin wind material beyond Keplerian rota- 
tion tend instead to lead to centrifugal mass ejection rather 
than a Keplerian disk. 

For much stronger fields, the region of rigid rotation be- 
comes more extended, and the MTD scenario can be viewed 



as becoming similar to the field-aligned rotation case (/3 = 0) 
of the RRM model developed here 9 . Although the disk ro- 
tation is rigid-body rather than Keplerian, the tendency for 
the bulk of the material to build up in the region near the 
Kepler radius means that the resulting line emission should 
develop a doubled-peaked profile that might be quite dif- 
ficult to distinguish from what is expected from a Keple- 
rian disk. Note, however, that such rigid-body disks seem 
unlikely to produce the long-term (years to decade) vio- 
let /red(V^R}_yjxiatiOTis often observed in Be-star emission 
lines jTelting et alJl994ll : such variations seem instead likely 
to be the result of long-t erm precession of elliptical o rbits 
within a Keplerian disk dSavoniie fc Heemskerklll993l) . As 
such, we do not believe that the RRM model is likely to be 
of general relevance to explaining Be star emission. However, 
as noted above, it does seem quite well-suited to explaining 
the rotationally modulated emission of Bp stars like a Ori E. 

On a concluding note, we draw attention to the fact 
that X-r ay flaring has been det ected in a Ori E by 
ROSAT JOroote fc Schmittl Eool) and s ubsequently by 
XMM-Newton iSanz-Forcada et all 120041) . iMullanl \200$) 
has argued that the flares originate from the B2 star itself, 
rather than from an unseen low-mass companion. If this is 
indeed the case, then we suggest a likely mechanism for the 
flare generation is thermal heating arising from magnetic re- 
connection. As we discuss in the Appendix, we expect the 
outer parts of the accumulation surface to undergo relatively 
frequent breakout events, during which stressed magnetic 
field lines will reconnect and release significant quantities of 
energy. We intend to exp lore this hypothesis f urther in a 
future paper (and see also lud-Doula et ai]|2004l) . 



8 SUMMARY 

We have presented a new Rigidly Rotating Magnetosphere 
model for the circumstellar plasma distributed around mag- 
netic early-type stars. By assuming that field lines remain 
completely rigid, and co-rotate with the star, we are able to 
find regions in the circumstellar environment where plasma 
can accumulate under hydrostatic equilibrium; in the gen- 
eral case of a tilted dipole field, these regions take the form 
of a geometrically-thin, warped disk, whose mean surface 
normal lies between the misaligned magnetic and rotation 
axes. When coupled with a quantitative description of the 
accumulation process, our treatment allows us to evaluate 
the density throughout the circumstellar environment, and 
thereby calculate observables such as emission-line spectra. 

This RRM model shows promise; even without a fine 
tuning of parameters, it reproduces the principal features 
of a Ori E, the archetype of the variable-emission helium- 
strong stars. In a forthcoming paper, we will investigate the 
extent to which the model can reproduce the more detailed 
aspects of this star - in particular, the strength of the emis- 
sion lines, as well as their shape, and the eclipse-line varia- 
tions seen in photometric indices. We also plan to examine 

9 Note, however, that while the MTD analysis emphasized the 
torquing role of the magnetic field, in the RRM model the rigid 
field also plays a crucial role in holding material down against 
a net centrifugal force that, for radii beyond the Kepler radius, 
exceeds the inward force of gravity. 
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whether the model can be applied to other magnetic early- 
type stars. 
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APPENDIX A: BREAKOUT OF 
ACCUMULATED MATERIAL 

Al Breakout Time 

In deriving the form of the accumulation surfaces, we have 
assumed an arbitrarily strong, rigid field. But in practice 
the finite magnitude of any stellar magnetic field means that 
there is a limit to the mass that can be contained against the 
centrifugal force; when the density gets too high, the mate- 
rial should break out from the field containment. Just prior 
to such an episode, the over-stressed magnetic field becomes 
distorted from its equilibrium configuration, sagging radially 
outward as it passes through the dense material in accumu- 
lation surfaces. Under these circumstances, the inward force 
arising from the tension in the distorted field lines barely bal- 
ances the net outward gravito-centrifugal force; therefore, we 
can construct an approximate condition for the occurrence 
of breakout by equating these two forces. Focusing our anal- 
ysis in this appendix on the simple case of an aligned dipole 
field (/3 = 0°), the breakout condition may be expressed as 



Pb 



GM, 



B 2 



47rh n 



(Al) 



where pb represents a breakout value for the peak density 
at radius r within the equatorial plane, and the scale height 
hm appears as the typical curvature radius of the distorted 
magnetic field lines, whose tension B 2 /iir generates the bal- 
ancing inward force. In analogy with eqn. 13011 we can define 
an associated breakout surface density <7b = v^Pb^m, where 
we have taken fj, m = 1 as appropriate to the equatorial plane. 
Using the Kepler radius rx (eqn. 1121 to scale both the local 
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radius (£ = r/rn) and the stellar radius (£* = R^/rn), we 
have for the usual dipole field scaling B ~ r -3 , 



or, solving for the breakout surface density, 



4ng t £ 4 (£ 3 - 1) 



(A2) 



(A3) 



Scaled in terms of typical parameters for a rotating, mag- 
netic B-star, we find a characteristic surface density for 
breakout 



8 gem 



.94 



(A4) 



where B 3 = B/10 3 G and g 4 = g*/W 4 cms -2 . 

For comparison, note that for a dipole field the surface 
density accumulation rate in eqn. (1 - > 1 has the scaling 



O-m(0 = 7? ■ 



£ 3 2tt7? 2 



(A5) 



where once more we have assumed fj, m = 1. Then for each 
scaled radius £ we can define a characteristic breakout time, 
tb(£) = (Jb/o"m- Casting the stellar gravity in terms of the 
surface escape speed and free-fall time, g* = v csc /2ts, the 
ratio of breakout to free-fall time becomes 



(A6) 



Here we have collected dimensional quantities in terms of 
a single, dimensionless magnetic confinement parameter for 
the accumulation surface 10 , 



BlRl ,„ e S|i?? 2 



M« 



1.6 x 10° 



M_io Vs 



(A7) 



with the latter equality giving a characteristic value in 
terms of scaled parameters R\2 = -R*/10 12 cm, M_io = 
M/10" 10 M yr~\ and v B = v CBC /10 8 cms" 1 . Noting that 
the free fall time tg = v CBC /2gt, = 10 4 s {v%/2g 4 ), the break- 
out time evaluates to 



t b (£) w 250yr 



M_ 10S4 £(£ 3 - i) 



(A8) 



where we have taken ypK I \i* ~ 2. As a typical example, 
corresponding roughly to values appropriate to a Ori E, let 
us take M- 10 = g 4 = 1, B 3 = 10, and R 12 = f, = 1/2, 
yielding a value 12.5 for the ratio factor in eqn. 1A8I . At 
a location equal to twice the Kepler radius, we then find a 
typical breakout time of it, (2) « 220 yr. 



A2 Mass in Accumulation Surface 

Let us next estimate the total mass in the equatorial accu- 
mulation surface after some elapsed time t since it was last 



emptied. The mass contained between inner radius n = £irK 
and outer radius r = £ fK is given by the integral 

m(t) w27rrfct^ ow(£K d £ 

PsMt/^Q.-i-V (A9) 

where the latter equality uses eqn. HA5II . Approximating the 
inner radius by the Kepler radius, £; « 1, we then find 



m(t) RiMtp^, 



(A10) 



Co ? 2 (eo 2 +Co + l) ' 

The latter equality assumes the outer radius is limited by 



breakout, and uses eqn. (IA6II to eliminate the explicit ap- 
pearance oft in terms of the time- variable outer radius £, {i). 
Over a long time the outer radius approaches the inner (Ke- 
pler) radius £ — * 1, with the total asymptotic disk mass 
approaching 



Mta^gV^ BlRlil^ 



(AH) 



3 6 S . ' 

where the latter equality uses the definition l|A70 to elim- 
inate the confinement parameter, mass loss rate, and ter- 
minal speed, and we have also eliminated the escape speed 
v csc in favor of the surface gravity g,. In terms of scaled 
parameters, this evaluates to 



1.5 x 1O" 8 M 



94 



(A12) 



Again adopting the above typical parameters for a Ori E 
- B3 = 10, £* = R12 = 1/2 and g 4 = 1 - we find moo as 
9.3 x 1O _8 M . 

If instead we consider a time when the outer radius 
happens to be at twice the Kepler radius, £ = 2, then by 
eqn. JA10I the total mass is reduced by an extra factor of 
3/28 = 0.107. For a Ori E, this now gives a total mass m w 
9.9 x 10 -9 M . As noted above, the associated breakout time 
for this twice-Kepler outer radius is about it, (2) ~ 220 yr. 

Overall, the picture from this analysis is that the outer 
parts of the accumulation surface should be subject to rela- 
tively frequent breakout events that empty mass from those 
regions. Over a longer time, rarer breakouts can occur from 
closer in, eventually even quite near the Kepler radius. This 
simple analysis formally assigns an arbitrarily long build-up 
time, and thus arbitrarily large mass build-up, to the Ke- 
pler ra dius itself. However, based on MHD simulations done 
so far JOwocki fc ud-Doulall2003D . it seems more likely that 
breakouts sufficiently close to the Kepler radius, i.e. with 
£ < 2, should be associated with a broader disruption of the 
overall field structure. This can lead to an emptying of mass 
throughout the accumulation surface, including the region 
around the Kepler radius itself. Following such a global evac- 
uation, the relative distribution of material at any given time 
is proportional to the wind accumulation rate, as detailed 
in Sec. El 



Note that this is closely related to the ' wind magnet ic confine- 
ment parameter' defined bv lud-Doula &; Owockil |2002|) . differing 
only by the order-unity substitutions Voo — > ti csc and B e q — > B„, 
where is the wind terminal speed, and B eq is the stellar sur- 
face field at the magnetic equator. 



